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1.0  INTRODUCTION 


In  the  last  three  decades,  the  use  of  devices  that  emit  radio  frequency  (RF)  fields 
has  increased  dramatically.  The  proliferation  of  RF  devices  has  been  accompanied  by 
increased  concern  about  ensuring  the  safety  of  their  use.  Throughout  the  world  many 
organizations,  both  government  and  non-government,  have  established  RF  safety 
standards  or  guidelines  for  exposure.  Because  of  different  safety  philosophies,  the  former 
USSR  and  some  of  the  Eastern  European  countries  have  more  stringent  safety  standards 
than  most  Western  countries.  The  maximum  permissible  exposure  in  some  Eastern 
European  and  Western  countries  differ  by  over  two  orders  of  magnitude.  These 
differences  have  raised  concerns  about  the  lack  of  uniformity  and  have  led  to  public 
concern  and  distrust  about  electromagnetic  fields  (EMF)  exposures  in  controlled  and 
uncontrolled  environments^  One  concern  about  these  differences  is  interpretation  of  data 
from  animal  studies.  Can  the  effects  observed  in  animals  exposed  to  EMFs  be  used  to  set 
human  exposure  standards?  Which  effects  are  important  for  human  safety  and  which  are 
not?  In  applying  safety  standards  the  issue  of  how  to  measure  EMF  absorption  in  the 
body  is  not  easily  solved.  Typically,  what  is  specified  in  the  safety  standard  is  field 
strength  or  power  density  that  should  not  be  exceeded.  These  limits  are  designed  to 
prevent  excessively  high  specific  absorption  rates  (SAR)  in  the  human  body.  Thus,  better 
methods  are  needed  to  properly  measure  EMF  absorption  and,  extrapolate  or  relate 
effects  observed  in  animals  to  those  expected  to  be  found  in  people.  The  resulting  data 
could  lead  to  modification  of  existing  safety  standards  or  estabUshment  of  new  safety 
standards. 

Accurate  dosimetry  represents  an  essential  element  of  the  research  in  determining 
the  biological  effects  of  EMF.  Electromagnetic  energy  is  absorbed  non-uniformly  in 
biological  tissues^^.  A  large  number  of  factors  such  as  a  body’s  shape  and  position,  as 
well  as  its  orientation  in  the  field  influence  these  non-uniform  distributions^^.  In  short, 
there  is  no  single  answer  to  die  question,  “How  much  electromagnetic  field  energy  will 
be  absorbed?”.  In  order  to  ensure  the  safe  use  of  EMF  emitting  devices,  a  number  of 
techniques  for  measuring  EMF  exposure  have  been  developed,  however,  they  all  have 
limitations.  Since  the  energy  absorbed  is  directly  related  to  the  internal  EMF  (that  is,  the 
EMF  inside  the  object,  not  the  EMF  incident  upon  the  object),  dosimetry  is  also 
interpreted  as  the  determination  of  internal  EM  fields.  First,  we  need  dosimetry  to 
determine  which  internal  fields  in  animals  cause  a  given  biological  effect.  Then  we  need 
dosimetry  to  determine  which  incident  fields  would  produce  similar  internal  fields  in 
people,  and  therefore  a  similar  biological  effect.  Dosimetry  is  needed  whether  the  effects 
are  produced  by  low-level  internal  fields  or  by  more  intense  fields  that  cause  body 
temperature  to  rise.  Since  internal  and  incident  EMF  can  be  quite  different,  depending  on 
the  size  and  shape  of  the  object,  its  electrical  properties,  its  orientation  with  respect  to  the 
incident  EMF,  and  the  frequency  of  the  incident  fields^;  the  development  of  theoretical 
dosimetry  modeling  techniques  and  powerful  computer  hardware  has  resulted  in 
dosimetry  modeling  as  a  principle  tool  in  determining  EMF  exposure. 

The  results  of  the  output  from  any  model,  including  Finite-Difference  Time- 
Domain  (FDTD),  are  questionable  if  the  model  has  not  been  validated  with  empirical 
data.  Validation  is  the  key  issue  in  earning  the  credibility  of  advanced  theoretical 
modeling  and  offers  a  unique  bridge  between  theoretical  and  experimental  dosimetry. 
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Since  applied  methods  in  experimental  dosimetry  require  the  implantation  of  thermal  or 
E-field  probes  mto  de^  body  organs,  it  is  not  feasible  to  perform  such  studies  in  human 
subjects,  therefore  vali^tion  attempts  of  the  computer  models  must  be  accomplished 
using  phantoms  (simplistic  geometry  and  tissue  equivalent  materials)  and  laboratory 
anmials.  Good  correlation  between  the  predicted  and  empirical  phantom  and  animal  data 
would  provide  increased  confidence  in  the  results  predicted  by  the  man  model^  A 

phantom  composed  of  a  sphere  is  a  usefol  and  convenient  method  to  validate 
FDTD  predictions  of  the  SAR  values^’^l  In  addition,  the  Mie  solution'^  which  provides 
m  exact  value  for  spheres,  can  be  used  to  measure  error  in  both  empirical  data  and  FDTD 
^ta.  Many  sophisticated,  high-resolution  digital  anatomical  models  that  are  being 
developed  call  for  validation  that  employs  not  only  simplistic  geometry  models,  but  also 
fee  actual  type  of  organisms  used  to  develop  these  models.  The  brain  of  fee  rat  offers  an 
Ideal  organ  for  companson,  due  to  fee  precision  of  stereotaxic  probe  placement,  as  well 
as  fee  ability  to  e^ily  confirm  probe  placements  in  tissue  sections  following 
experimentation  , 

The  facilities  at  fee  Radiofi-equency  Radiation  Branch  (HEDR),  Air  Force 

Research  Laboratory,  for  accomplishing  these  empirical  measurements  are  described 
below. 


2.0  EXPERIMENT  AT.  DOSIMETRY 

Measurements  of  fee  external  fields  for  fee  evaluation  of  exposure  conditions  are 
usually  considered  a  part  of  expaimental  dosimetry.  To  place  the  measurements  of 
external  fields  vs.  internal  fields  in  proper  perspective,  one  must  realize  that  fee 
instruments  and  methods  for  measurements  of  internal  fields  are  not  suitable  for  use  on 
humans  except  in  hyperthermic  treatment  for  cancer  feerapy^^  Therefore  to  prevent 
potentially  hazardous  exposures  to  humans,  fee  external  fields  are  measured.  Data  on  fee 
intensities  of  external  fields  feat  potentially  constitute  a  hazard  to  man  are  obtained 
toou^  experimental  evaluation  of  fee  intensities  of  the  internal  fields  and  fee  resulting 
determination  of  external  fields  which  produce  similar  internal  fields 

(SAR)  m  man. 

2.1  Exposure  Facilities 

During  many  radio  fi-equency  (RF)  exposures,  it  is  important  to  know  the  incident 

power  density.  There  are  many  means  of  determining  this  power  density.  The  three  most 

common  used  are  a  RF  survey  meter  in  conjunction  wife  a  compatible  probe,  a  calibrated 
power  sensor  and  an  antenna  wife  a  well  defined  capture  area,  and  temperature 
fefferences  (AT)  measurements  using  an  infi-ared  camera  and  carbon-loaded  Teflon  (See 
Tables  1-3  and  Fig.  1-5).  For  fee  fi-equency  range  of  30  kHz  to  40  GHz,  the  favorite 
method  is  using  the  RF  survey  equipment. 

There  are  a  few  mmufacturers  of  RF  survey  equipment.  Narda  Microwave 
manufactures  fee  primary  instrument  used  to  measure  electromagnetic  fields.  The  two 
meters  used  are  fee  8617  malog  meter  wife  appropriate  probe  or  fee  8718  digital  meter. 
Usmg  either  meter  is  relatively  easy.  The  probe  and  meter  are  calibrated  together  and 
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zeroed  in  a  zero  field  environment.  The  probe  is  then  placed  in  the  RF  field  to  be 
measured.  The  probes  are  isotropic,  but  the  best  practice  is  to  face  the  probe  head  toward 
the  RF  source.  There  are  three  elements  in  the  probe  head,  so  in  order  to  best  align  the 
elements,  the  probe  is  rotated  through  360°  and  measurements  recorded  every  45°.  These 
numbers  are  then  averaged  and  the  standard  deviation  is  calculated.  Next,  add  the 
standard  deviation  to  the  average  and  apply  the  correction  factor  for  the  probe.  Using  a 
spreadsheet  makes  the  calculations  very  quick  and  easy.  For  lower  and  higher  fi-equency 
fields,  other  techniques  were  developed. 

Millimeter  waves  are  not  easily  measured  using  the  RF  survey  equipment.  A 
method  using  a  power  sensor  and  a  capture  device  was  developed.  A  combination  of  an 
HP/Agilent  W8486A  waveguide  power  sensor  and  Millitech  standard  gain  pyramidal 
horn  connected  to  an  Agilent  437B  power  meter  is  used  to  measure  file  fields  fi'om  75  - 
110  GHz.  The  USAF  Precision  Measurements  Equipment  Laboratory  (PMEL)  calibrates 
both  the  power  sensor  and  power  meter.  The  effective  area  (Ae)  of  the  horn  is  determined 
fi’om  the  actual  gain  of  the  hom.  For  very  small  spot  sizes,  an  open-ended  waveguide 
section  can  be  used  in  lieu  of  a  hom.  The  Ae  however,  is  not  easy  to  determine.  The  only 
reliable  method  of  determining  the  correction  factor  is  to  do  a  side-by-side  comparison 
with  a  known  measurement  device.  The  open-ended  waveguide  is  very  fi-equency  and 
angle  sensitive.  Attenuators  can  be  placed  in  line  after  the  hom  to  extend  the  useful 
power  measurement  range.  For  very  high  power  densities,  other  methods  must  be  used. 

Numerous  materials  were  investigated  as  tools  to  over  come  the  restriction  of 
power  density  or  spot  size.  These  materials  were  looked  at  with  an  infi-ared  camera  to 
measure  temperature  rise  verse  power  density.  RF  cloth  or  carbon-impregnated  cloth  was 
first  tried.  This  material  responds  to  RF  exposure  with  a  very  uniform  temperature  rise, 
but  with  no  significant  thermal  mass,  it  was  very  susceptible  to  environmental  conditions. 
Carbon  loaded  Teflon  was  tried  next.  This  material  has  very  uniform  characteristics  and 
has  enough  mass  to  be  minimally  effected  by  environmental  conditions.  Extensive 
temperature  rise  measurements  using  an  infi-ared  camera  were  performed  in  an 
environmentally  controlled  chamber.  Power  density  was  measured  using  the  waveguide 
sensor  described  earlier.  A  series  of  points  were  generated  and  an  equation  developed  to 
correlate  the  temperature  rise  with  power  density.  The  carbon-loaded  Teflon  is  also  very 
useful  in  mapping  fields  that  have  interference  patterns.  The  maximum  and  mitiimnm 
points  are  easily  determined.  From  the  interference  patterns,  the  average  delivered  power 
can  be  calculated.  Other  uses  for  carbon-loaded  Teflon  are  still  being  investigated.  More 
accurate  and  defined  calibration  procedures  are  being  developed.  With  time,  this 
technique  should  prove  very  useful. 

These  are  only  a  few  types  of  measurements  that  can  be  made.  Dipole  antennas 
and  H-field  loops  are  some  other  common  methods  that  have  been  used^^.  Other 
techniques  may  allow  us  to  make  more  accurate  measurements. 

All  power  measurement  equipment  and  blackbody  source  cahbration  is  performed 
by  or  certified  by  PMEL  located  at  Brooks  City-Base  under  the  31 1  Mission  Support 
Group/Logistics  (31 IMSG/LG). 
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Table  1.  RF  Measurement  Resources 


RESOURCE 

MAKE 

MODEL 

RANGE 

ACCURACY 

RF  Power  Meter 

Hewlett  Packard 

HP437B 

Determined  By 
Probe  Used 

Absolute  Mode  +/- 
.02dB  or  +/-.5  % 

RF  Power  Meter 

Narda 

8718 

Determined  By 
Probe  Used 

Determined  By 
Probe  Used 

RF  Power  Sensor 

Hewlett  Packard 

W8486A 

75  GHz  to  110 
GHz 

Power  Linearity 
±D 

RF  Power  Sensor 

Hewlett  Packard 

R8486A 

26.5  GHz  to  40 
GHz 

Power  Linearity: 
±3% 

RF  Power  Probe 

Narda 

8723 

300MHz  to  50GHz 

+L25/-3.0  dB  (0.3 
GHz  to  50  GHz) 

±  1.25  dB(l  GHz 
to  50  GHz) 

RF  Power  Probe 

Narda 

8762 

300KHztolGHz 

±  dB(3  MHz  to 
300  MHz)  ±2  dB 
(300  kHz  to  1 
GHz) 

Standard  Gain 
Horn 

Millitech 

SGH-IO-RPOOO 

75-llOGHz 

Gam=23.2  db  @ 
93.2  GHz 

Waveguide 

Attenuator 

Dorado 

FAlO-10 

75-llOGHz 

±1.5db 

IR  Camera 

Flir 

ThermaCAM® 

S60 

-40°C  -1500  °C 

±2°Cor±2% 

IR  Camera 

Flir 

ThermaCAM®  SC 
3000 

-2b°C  to  1500°C 

±  l°Cor±  1% 

IR  Camera 

Indigo 

MerUn  Mid 

0  to  350°C 

±2°Cor±2% 

IR  Camera 

Amber 

Radiance- 1 
Thermographic 

Determined  By 
Black  Body  Cal. 

Determined  By 
Black  Body  Cal. 

Open-ended 

Waveguide 

Fabricated  In- 
house 

N/A 

75-llOGHz 

Determined  By 
Comparison 

Blackbody  Source 

Mikron 

340 

-20  to  150‘'C 

±.3°C 

Blackbody  Source 

Mikron 

345X8 

-10.00  to  148.9°C 

<  20mk 

A 


Table  2.  Available  Exposure  Systems 


Freq 

(MHz) 

H 

Bldg 

Manufacture 

Comments 

cw 

Chi 

Amplifer  Research 

dipole  antennas 

cw 

HI 

Chi 

Southwest  Research 
Institute 

Solid  State 

3k 

1161 

Ch2 

Polarity  Inc 

Delivery 
cy  2004 

915 

75k 

1162 

1130 

Microdry 

2450 

Ik 

CW 

1162 

Cober 

Portable 

2450 

Ik 

CW 

1162 

MCL 

Portable 

35000 

75 

CW 

1162 

1133 

Applied 

Elecromagnetics 

Dual  Freq 
Cabinet 

95000 

50 

CW 

1162 

1133 

Applied 

Elecromagnetics 

Dual  Freq 
Cabinet 

95000 

100 

CW 

1125 

CPI 

UWB 

1162 

1077 

Various 

750-975 

■■ 

cw 

1183 

Ch9 

MRC-113 

Removed  from 
Van 

2070 

■E3I 

cw 

1183 

Ch8 

Cober 

8300 

cw 

1183 

Ch9 

Cober 

35000 

3500 

_ 

1183 

Ch7 

General  Electric 

Scheduled  for 
DRMO 

2070 

1185 

Ch3 

Cober 

Part  of  PPES 

2800 

2M 

1185 

Ch3 

Cober 

Part  of  PPES 

■a 

IM 

1185 

Ch3/ 

Ch4 

Cober 

Part  of PPES 

■essehi 

700 

700k 

1185 

Ch3 

Cober 
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Figure  1.  Applied  Elecromagnetics  35/94  GHz  System. 


Figure  2.  Peak  Power  Emission  System  (PPES) 


Envtronmentai 

Control 


22x32x22  10x20x10  >50  MHz  18-30°C±0.5°C 

12x20x10  10x18x9  >  1  GHz 

24  X  40  X  10  22  X  38  X  9  >1  GHz 


8x12x8  6x10x6  >  1  GHz  18-30°C±0.5°C 
10  X  14  X  8  8x12x6  >1  GHz  10-35°C  ±0.5®C 


20x24x10  20x24x10 


8x10x8  7x9x7 


>1  GHz 


10-35“C±0.5‘’C 


10-35'’C±0.5“C 


10-35“C±0.5“C 


15-30“C±1.0®C 


None 


15-30'’C±1.0“C 


2.2  Measuring  SAR 


Experimental  dosimetry  includes  the  following  areas;  measurements  of  total  SAR 
and  its  distribution  in  live  or  dead  experimental  animals,  measurements  of  average  SAR 
and  SAR  distribution  in  models  of  animals  and  humans  (phantoms),  and  development  of 
exposure  devices  for  animals  and  in  vitro  preparations  in  which  quantification  of 
exposure  conditions  and  internal  fields  is  facilitated’’.  Experimental  dosimetry  plays  a 
crucial  role  in  the  control  of  animal  exposures  and  quantification  of  the  data  fi'om  animal 
experimentation.  It  enables  comparison  of  experiments  conducted  in  various 
laboratories,  under  different  exposure  conditions,  and  for  different  species.  It  also  plays 
an  important  role  in  validation  of  the  calculations  from  theoretical  dosimetry.  The 
dosimetric  data  can  be  used  to  extrapolate  the  results  of  animal  experiments  to  human 
exposures  where  the  physical  scaling  principle  is  employed^.  Generally,  accepted 
methods  of  measurement  of  SAR  include  measurement  of  the  rate  of  temperature  rise 
within  the  exposed  object  or  measurement  of  the  internal  E-field  strength.  Temperature 
rise  may  be  characterized  by  a  whole-body-averaged  (calorimetric)  measurement,  a  point 
measurement  (via  a  thermometer  implanted  in  the  body  being  exposed),  or  infi-ared 
thermographic  camera  analysis  of  bisected  phantom  models  that  have  been  exposed  to 
large  RF  fields’*.  The  internal  electric  field  strength  may  be  measured  by  an  implantable 
E-field  probe’®. 

2.3  Temperature  Measurement  Techniques 

An  important  advantage  of  the  temperature  measurement  technique  is  that 
temperature  is  a  scalar  value,  making  it  easier  to  produce  probes  of  very  small  size  having 
short  time  constants^.  One  major  difficulty  is  the  need  for  high  power  RF  sources  in 
order  to  achieve  large  enough  SAR  values  for  accurate  measurement,  i.e.  <  1  W/kg. 
Another  important  drawback  is  that  the  time  derivative  can  only  be  determined  by 
temperature  measurements  at  two  discrete  times,  whereby  the  tiiermodynamic  state  prior 
to  exposure  must  be  known  and  should  preferably  be  in  equilibrium’®.  Since  the 
exposure  alters  thermal  equiUbrium,  means  must  be  provided  to  re-establish  equilibrium 
after  each  measurement  cycle,  which  is  usually  time  consuming. 

Average  whole  body  SAR  may  be  measured  using  calorimetric  methods.  Such 
methods  have  been  used  predominantly  with  small  animals  or  animal  models*’  but 
also  to  measure  SAR  in  a  full-size  human  model’’.  Although  the  calorimetric 
determination  of  the  energy  stored  in  the  model  once  inside  the  calorimeter  is  quite 
accurate  per  se,  the  overall  system  accuracy  in  terms  of  SAR  is  limited  by  how  closely 
the  test  object  models  the  actual  object  and  by  the  amount  of  irradiation-induced  heat 
that  escapes  fi-om  the  object  that  is  not  measured.  This  method,  moreover,  requires 
sufficient  time  for  the  thermal  equilibration  processes  to  occur  and  requires  sufficient 
energy  deposition  in  the  test  object  to  produce  a  calorimeter  output  signal  that  is  enough 
above  baseline  to  be  measurable. 

However,  the  Dewar-flask  method  of  calorimetry  is  a  relatively  straightforward 
way  of  determining  the  whole-body  average  SAR  of  small-bodied  animals”’  ’*.  The 
calorimetric  technique  of  determining  a  whole-body  average  temperature  requires  that 
the  cadaver  be  immersed  in  a  Dewar  flask  containing  a  medium,  such  as  water,  at  a 
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known  temperature;  then  the  temperature  of  the  carcass,  following  irradiation,  can  be 
determined  by  calculating  the  final  temperature  of  the  carcass. 

Non-perturbing  temperature  probes  (small,  high  resistance  thermistors  with  high 
resistance  leads)  can  be  used  successfully  to  make  SAR  measurements  especially  point 
measurements  in  very  small  structures  (e.g.  small  locations  in  animars  brain).  The 
minimum  requirements  are  that  the  temperature  sensor  and  associated  leads  should  be 
nonperturbing  to  the  EMF,  and  the  SAR  should  be  large  enough  to  produce  a  measurable 
temperature  rise  during  a  period  of  about  30  seconds.  The  first  requirement  can  be 
satisfied  by  using  highly  resistive  material  or  fiber  optics,  instead  of  metal  components, 
for  the  temperature-sensing  element  leads^^.  The  second  requirement  entails  the 
measurement  of  SARs  no  lower  than  a  few  W/kg.  This  lower  limit  exists  because  the 
resolution  of  most  temperature  probes  is  typically  0.01  to  0.1°C,  and  the  longest  practical 
duration  of  irradiation  that  allows  reasonably  accurate  SAR  measurement  is  typically 
between  20  and  30  seconds'^  because  longer  exposures  allow  heat  loss  to  occur. 

Our  instrument  of  choice  is  an  in-house  developed  temperature  measurement  and 
acquisition  system  (TMAS)(See  Fig.  6  and  Table  4).  It  was  designed  around  the  BSD 
non-perturbing  thermister  probe  developed  by  Ronald  R.  Bowman^^  In  addition  to  the 
BSD  probe,  the  system  consists  of  an  8-channel  interface  box^^  a  National  Instruments 
acquisition  card,  and  Lab  VIEW  based  software  run  on  a  notebook  computer. 

The  interface  consists  of  a  highly  stable  precision  constant  current  source 
providing  excitation  for  the  thermister.  The  voltage  developed  across  the  thermister  is 
fed  into  a  very  high  impedance  instrumentation  amplifier.  The  signal  is  scaled  to  match 
the  A-D  card  and  then  low-pass  filtered  at  10  Hz.  This  signal  is  then  cabled  to  the  A-D 
card  located  in  the  computer.  The  front  panel  provides  input  jacks  for  the  probes  and  a 
power  LED.  The  rear  panel  contains  the  On/OfF  switch  power  connector  and  an  output 
connector.  The  interface  is  powered  by  common  120VAC. 

The  Software  was  developed  using  the  Graphical  programming  language 
Lab  View.  It  quantatizes  the  interface  outputs  at  400  samples/sec,  averaging  80  samples  5 
times/sec.  This  provides  5  data  points/sec  max,  which  is  based  on  the  time  constant  of 
the  probe.  The  software  linearizes  the  thermister  output  using  the  Steinhart  and  Hart’° 
equation  for  thermisters.  Coefficients  for  the  equation  are  derived  during  the  calibration 
process.  The  software  provides  spreadsheet  compatible  data  that  includes  date,  time  and 
sample  number  from  start  for  each  data  point.  The  computer  display  provides  all 
operational  controls,  data/file/path  display  and  real-time  charting. 

A  thermographic  camera  can  be  used  to  measure  temperatures  and  ultimately 
SAR  across  the  visible  surface  of  an  object'^’ Since  the  camera  is  non-invasive  it  ran 
be  ^ed  in  addition  to  other  measurement  techniques.  The  method  involves  the  use  of  a 
radiometric  camera  for  recording  the  rate  of  rise  of  temperature  in  a  plane  that  bisects  an 
entire  model  (phantom)  or  a  cadaver  of  a  biological  subject  under  study.  The  temperature 
distribution  before,  and  immediately  after  brief,  high-power  irradiation,  is  observed  on 
the  precut  surface  on  each  half  of  the  bisected  model.  This  is  done  to  prevent  cooling  by 
evaporation  or  flow  of  the  wet  synthetic  tissue  out  of  its  shell  (usually  composed  of 
S3mthetic  fat  or  rigid  plastic  foam).  This  procedure  is  also  used  to  measure  heating 
patterns  on  the  surface  of  subjects  exposed  to  RF  fields^*.  Infrared  thermometry  offers  at 
least  three  advantages:  1)  temperature  measurements  are  recorded  in  real  time,  by  a 
system  with  fast  effective  response  time,  thus  eliminating  any  temporal  error  due  to  the 
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finite  response  time  of  thermistors  or  thermocouples,  2)  the  sensor  array  of  the  IR  camera 
allows  the  simultaneous  sampling  of  multiple  discrete  areas,  and  (3)  no  sensor  is  placed 
in  contact  with  the  tissue,  which  reduces  the  chance  of  artifact.  This  latter  advantage 
assures  that  measurements  are  obtained  from  the  site  that  reaches  the  maximal 
temperature,  thus  eliminating  error  due  to  less-than-optimal  placement  of  a  single  probe. 


Figure  6.  Tenqjerature  measurement  and  acquisition  system  (TMAS). 


Table  4.  HEDR  TMAS  Specifications 


Accuracy: 

0.1  deg  C 

Resolution: 

O.Oldeg  C 

Stability: 

+/-.  05deg  C,  Tested  to  one  month 

Channels: 

8 

Range: 

18-60degC 

Software: 

Lab  VIEW  Custom  Application 

Power: 

^i20VAC 
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2.4  Electric  Field  Measurement  Techniques 

Since  thermal  as  well  as  non-thermal  interactions  of  RF  fields  with  organisms  are 
related  to  the  internal  E-field  intensity,  E- field  sensors  have  been  developed^^.  In  fact,  the 
simplest  way  of  measuring  electric  field  is  to  use  a  monopole  or  dipole  anteima  at  the  end 
of  coaxial  line  that  is  connected  to  a  detection  and  measuring  system. 

Miniature  isotropic,  three-dimensional  diode-loaded  E-field  probes  with  high 
impedance  feed  lines,  which  have  been  commercially  available  for  a  number  of  years, 
have  been  used  to  measure  SAR  distributions  in  phantom  models  and  in  living, 
anesthetized  animals^®.  This  type  of  probe  can  be  inserted  into  a  sample  in  the  same 
manner  as  the  non-perturbing  temperature  probes.  These  probes  measure  the  intensity  of 
the  internal  E-field  at  the  location  of  the  probe. 

2.5  Power  Difference  Method 

Whole-body  dosimetry  may  be  performed  using  a  Power  Difference  method  that 

is  limited  to  closed  exposure  systems  such  as  waveguides,  TEM-mode  transmission  lines, 
cavities,  etc.  In  this  method,  power  meters  are  attached  to  the  exposure  systems  by 
directional  couplers  to  measure  the  incident,  reflected,  and  transmitted  power.  The  power 
absorbed  by  the  biological  target  (whole-body  absorption  rate)  is  then  measured  by 
subtracting  the  values  of  the  reflected  and  transmitted  power  from  the  incident  power. 

For  exposure  systems  using  live  ammals,  the  values  of  these  three  measurements  vary  as 
the  animals  change  position  and  posture  because  absorption  is  a  function  of  the  animal’s 
orientation  relative  to  the  field  polarizations.  It  is  thus  desirable  to  have  a  data  collection 
system  that  can  integrate  the  measured  power  difference  and  determine  a  temporal 
average,  usually  for  the  duration  of  exposure.^  ^  The  advantage  of  this  method  is  that  it 

is  capable  of  providing  on-line  instantaneous  measurements  of  the  whole-body  absorption 
rate  in  live  animals.  As  such,  the  method  is  useful  for  continuously  monitoring 
absorption  changes  created  by  alterations  in  an  animal’s  position,  posture,  etc.  There  are, 
however,  some  problems.  In  cases  involving  a  multiple-animal  exposure  within  the  gamp 
closed  system,  the  measured  absorption  rate  represents  the  sum  total  for  all  animals. 

How  this  energy  is  distributed  among  the  animals  within  the  system  caimot  be 
determined.  The  method  is,  however,  optimal  for  single-animal  exposures. 


3.0  SAR  MODELING 


The  last  decade  has  witnessed  significant  progress  in  numerical  dosimetry  for 
human  exposure  to  RF  fields.  With  the  development  of  the  FDTD  method,  it  is  possible 
to  predict  electric  field  strength  and  SAR  values  in  various  biological  tissues  under  a 
wide  range  of  exposure  conditions  in  phantoms  and  laboratory  animals,  as  well  as 
humans.  Permittivity  values  play  a  dominant  role  in  the  overall  consideration  of 
interaction  between  EMF  and  matter  and  in  related  applications  including  EM  dosimetry. 

The  rigorous  analysis  of  a  realistically  shaped  heterogeneous  model  for  humans 
or  experimental  animals  is  a  difficult  theoretical  task.  Because  of  the  difficulty  of 
solving  Maxwell’s  equations,  which  form  the  basis  of  this  analysis,  a  variety  of  special 
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models  and  techmques  have  been  used,  each  valid  only  in  a  hmited  range  of  frequencies 
or  other  parameters. 

With  the  availability  of  improved  computer  technology,  various  numerical 
techmques  have  been  developed  for  the  solution  of  EMF  interactions  with  objects  having 
irregular  geometries  and  heterogeneous  electrical  composition.  Due  to  mathematical 
complexities  involved  in  calculating  SAR,  a  combination  of  techniques  has  been  used  to 
obtain  SAR  for  various  models  as  functions  of  frequency.  Each  of  these  techniques 
provides  information  over  a  limited  range  of  parameters.  Combining  the  information 
thus  obtained,  gives  a  reasonably  good  description  of  SAR  as  a  function  of  frequency 
over  a  wide  range  of  frequencies  for  a  number  of  useful  models. 

From  a  wide  array  of  methods,  including  the  extended-boundary  condition 
method  (EBCM)  and  iterative  extended-boundary  condition  method  (lEBCM),  method  of 
moments  (MOM),  finite  element  method  (FEM),  finite-element  time-domain  method 
(FETD),  generalized  multipole  technique  (GMTX  the  conjugate  gradient-fast  Fourier 
transform  (CG-FFT),  volume-surface  integral  equation  method  (VSIE),  quasi-static 
admittance  and  impedance  methods,  the  FDTD  method^^'^’  has  become  the  most  widely 
used  for  bioelectromagnetic  applications  in  the  range  of  a  few  MHz  to  several  GHz.  An 
extension  of  the  FDTD  method,  the  frequency-dependent  finite-difference  time-domain 
method  ((FD)^TD)^*’  enables  broadband  bioelectromagnetic  simulations  by  including 
the  effect  of  the  frequency  dispersion  of  the  tissues.  This  method  has  been  used  to 
calculate  SAR  and  cxirrent  distributions  in  the  body  resulting  from  ultra-short  plane  wave 
pulses  with  bandwidths  of  the  order  of  1  GHz. 


3.1  Finite-Difference  Time-Domain  Code 

Kunz  and  Luebbers^^  originally  developed  an  FDTD  program  to  predict  SAR 
distributions  in  different  digital  models  (sphere,  rat,  and  human).  This  FDTD  code  was 
modified  and  digital  anatomical  models  were  developed  jointly  by  U.S.  Naval  Health 
Research  Center  Detachment  and  U.S.  Air  Force  Research  Laboratory,  Brooks  City-Base, 
Texas.  The  permittivity  properties  of  each  of  the  tissue  types  are  set  according  to  data  and 
approximations  published  by  Gabriel'^®.  The  modified  code  reads  the  anatomical  model 
files  and  outputs  a  number  of  files  that  include  3-D  normalized  SAR;  mean,  minimuTn  and 
maximum  SARs  for  each  tissue  type;  and  each  Z-plane  slice.  Finally,  there  is  an  extensive 
log  file  and  all  file  names  are  automatically  constructed  to  reflect  run  parameters.  The  code 
has  been  parallelized  using  the  message-passing  interface  (MPI)  library,  which  allows  for 
larger  and  more  complex  data  sets  to  be  modeled.  The  advantage  of  using  the  MPI  is  that 
the  code  can  run  on  parallel  computer  systems  composed  of  networks  of  computers.  In 
order  to  run  a  human  model  with  1  mm  resolution,  we  had  a  Linux-based  Beowulf  parallel 
processing  system  with  100  processors  and  a  total  of  104  GB  of  RAM  developed  (Fig.  7). 

The  FDTD  algorithm  was  first  established  by  Yee'*^  as  a  three  dimensional 
solution  of  Maxwell's  curl  equations; 


Vx£  =  -— =  -//— 
dt  ^  dt 


VxH  =  aE  +  £ 


dt 
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In  principle,  a  volume  of  space  containing  any  object  or  collection  of  objects  is 
subjected  to  an  electromagnetic  disturbance,  FDTD  then  solves  for  the  fields  throughout 
the  volume  as  a  function  of  time.  The  FDTD  method  was  later  improved^^"^^. 

It  uses  a  geometry  mesh,  usually  of  rectangular  box-shaped  cells  (voxels),  which 
is  readily  developed  from  CT  or  MRI  scans  of  real  human  beings  or  animals  The 
constitutive  parameters  for  each  cell  may  be  set  independently  so  that  objects  having 
irregular  geometries  and  heterogeneous  electrical  composition  can  be  analyzed.  The 
FDTD  method  has  been  used  for  myriad  applications  including  calculating  SARs  and 
induced  currents  in  the  human  body  for  plane-wave  exposures,  exposure  to  leakage  fields 
of  parallel-plate  dielectric  heaters,  exposure  to  pulse  EMF,  annular  phased  arrays  of 
aperture  and  dipole  antennas  for  hyperthermia,  coupling  of  cellular  telephones  to  the 
head,  exposure  to  RF  magnetic  fields  in  magnetic  resonance  imaging  (MRI)  machines 
and  exposure  to  power  line  fields. 

The  FDTD  method  uses  simple  central-difference  approximations  to  evaluate  the 
space  and  time  derivatives  using  a  time  stepping  procedure.  Inputs  are  time-sampled 
analog  signals.  The  region  being  modeled  is  represented  by  two  interleaved  grids  of 
discrete  points.  One  grid  contains  the  points  at  which  the  magnetic  field  is  evaluated. 

The  second  gnd  contains  the  points  at  which  the  electric  field  is  evaluated. 

Values  of  permeability  (p),  relative  dielectric  (s’)  and  effective  conductivity  (ct) 
assigned  to  each  field  component  in  each  cell  define  the  position  and  electrical  properties 
of  the  scatterer.  These  parameters  can  have  different  values  for  different  field 
orientations  permitting  amsotropic  materials  to  be  modeled.  Their  values  can  also  be 
adjusted  at  each  time-step  depending  on  conditions  making  it  easy  to  model  nonlinear 
materials.  Because  the  basic  elements  are  cubes,  curved  surfaces  on  a  scatterer  must  be 
"'stair  casQd'\  For  many  configurations  this  does  not  present  a  problem.  However  for 
configurations  with  sharp,  acute  edges,  an  adequately  stair  cased  approximation  may 
require  a  very  small  grid  size.  This  can  significantly  increase  the  computational  size  of 
the  problem.  Surface  conforming  FDTD  techniques  with  non-rectangular  elements  have 
been  introduced  to  combat  this  problem.  One  of  the  more  promising  of  these  techniques, 
which  permits  each  element  in  Ae  grid  to  have  an  arbitrary  shape,  is  referred  to  as  the 
Finite  Volume  Time  Domain  (FVTD)  method‘s. 

As  more  powerful  computers  became  widely  available,  use  of  the  FDTD  method 
has  increased  exponentially^^.  Its  use  in  bioelectromagnetics  has  been  increased  as 
well"**'^®.  Unlike  the  alternatives,  the  method  of  moments  (MOM)^’  and  finite  element 
method  (FEM)  ,  die  FDTD  method  is  scalable,  i.e.,  the  CPU  time  behaves  linearly  with 
problem  size  N.  The  MOM  and  FEM  methods  require  matrix  inversions  (albeit  sparse 
matrices  in  FEM)  and  thus  scale  as  N^  The  MOM  method  is  primarily  useful  for 
problems  with  conducting  surfaces,  but  is  difficult  to  apply  to  permittivity  problems  of 
interest  here.  Its  practical  limitation  is  to  systems  with  less  than  10^  cells^^.  The  FEM 
method,  with  its  irregular  cell  stracture,  is  difficult  to  parallelize  efficiently.  The  FDTD 
method  with  its  rectangular  cell  structure  is  easy  to  parallelize,  and,  in  the  case  of  our 
problem,  is  compatible  with  the  cellular  data  formats  of  the  animal  and  human  anatomical 
models.  The  main  disadvantages  of  FDTD  are  object  resolution  and  absorbing  boundary 
conditions  (also  true  in  FEM),  but  s^histicated  versions  of  the  FDTD  method  have  been 
developed  to  handle  these  problems^. 
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We  also  have  the  XFDTD  version  6.0  Bio-Pro  from  Remcom  running  on  dual 
Pentium  HI  (1  GHz)  processors  with  4  GB  of  RAM. 


Figure  7.  Linux-based  Beowulf  parallel  processing  system. 


3.2  Models 

Over  the  years,  experimental  phantoms  have  been  developed  to  understand  coupling 
of  EM  fields  to  models  of  the  biological  systems^^.  Early  analyses  were  based  on  spheres, 
cylinders,  spheroids,  ellipsoids,  and  large-block  models  (cubical  mathematical  cells 
arranged  in  file  shape  of  a  human  or  animal  body).  Although  these  models  were  relatively 
crude  representations  of  the  size  and  shape  of  the  human  body,  experimental  results  show 
that  calculations  of  the  average  SAR  agree  reasonably  well  with  empirical  values^. 
Calculations  of  the  local  distribution  of  the  SAR,  even  though  much  more  difficult,  are  now 
becoming  possible.  Simple  homogeneous  models  have  been  used’^.  Though  these 
simplistic  models  are  incapable  of  giving  accurate  SAR  distributions,  they  do  allow  use  of 
actual  near-field  sources  such  as  wireless  telephones.  And,  they  have  been  shown  to  give 
peak  1-g  SARs  needed  for  SAR  compliance  testing  within  an  uncertainty  of  ±  20 

Concomitant  with  the  development  of  homogeneous  models,  the  resolution  of 
heterogeneous  anatomically  based  models  of  the  human  body  has  been  increasing.  Initially 
based  on  anatomical  sectional  diagrams,  more  recent  versions  are  based  on  MRI  scans  of 
living  humans^^'^^  or  cadavers  such  as  "visible  man"  and  "visible  woman"  developed  by 
National  Library  of  Medicine^  *(See  Fig.  9).  Development  of  models  suitable  for  FDTD 
dosimetric  calculations,  while  straightforward,  is  not  trivial.  MRI  and  CT  scans  provide 
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voxel  maps  of  density,  but  many  of  the  tissues  have  the  same  or  similar  densities,  and  many 
of  the  regions  outside  of  the  major  organs  require  a  detailed  understanding  of  anatomy  to 
determine  what  tissues  are  present,  e.g.,  fat,  muscle,  bone.  Even  when  some  automatic 
tissue  definition  can  be  used,  the  bulk  of  the  task  falls  on  manual  tissue  recognition  by  a 
trained  anatomist.  Moreover,  since  FDTD  calculations  use  the  tissue  electrical  properties, 
i.e.,  dielectric  value  and  conductivity,  and  mass  density  at  each  voxel  location,  knowledge 
of  these  properties  for  each  cell  is  also  required.  With  the  advent  of  larger  memory  in  PCs, 
computer  workstations,  and  parallel  processors,  it  is  possible  to  run  anatomically  based 
whole-body  models  with  resolutions  on  the  order  of  a  millimeter  and  even  sub-millimeter 
resolution  models  for  partial  body  exposures  (e.g.  head  and  neck)  for  exposures  such  as 
wireless  telephones  .  Anatomically  based  models  of  the  common  laboratory  animals  such 
as  the  rat,  goat,  and  rhesus  monkey  have  also  been  developed^^  and  used  for  dosimetric 
calculations  because  of  the  importance  of  such  subjects  for  laboratory  experiments^*^’ 
Development  of  the  anatomical  models  of  the  Sprague-Dawley  rat  and  man  has 
been  described  in  details'*’  Briefly,  magnetic  resonance  imaging  (MRI)  was  used  to 
acquire  axial  scans  of  the  rat.  Images  of  the  man  were  obtained  fi’om  the  Visible  Human 
project  (National  Library  of  Medicine).  Initial  computer-segmentation  of  tissue  types  on 
the  man  images  was  completed  by  CieMed  (a  collaboration  between  National  University  of 
Singapore  and  Johns  Hopkins  University).  The  original  x,  y,  z  voxel  dimensions  were  .39 
X.39  X  3.0  mm  for  the  rat  and  1 .0  x  1 .0  x  1 .0  mm  for  the  man.  Each  voxel  was  color-coded 
and  assigned  a  tissue  type.  In  addition  to  these  models,  we  have  produced  computer 
anatomical  models  of  a  rhesus  monkey,  phantom  monkey,  pigmy  goat,  and  rat  (Fig.  8). 
These  models  are  available  via  anonymous  FTP  fi-om  starview.brooks.af.mil/EMF. 


Figure  8.  Models  of  A)  Rat,  B)  Goat,  C  &  D)  Monkey,  and  E)  Phantom  Monkey. 
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Figure  9.  Visible  Man  Model 


We  also  are  running  Varipose  from  Remcom  which  allows  us  to  reposition  the 
visible  human  mesh.  This  software  was  developed  through  the  AFRL-sponsored  Small 
Business  Innovative  Research  program. 

Prof.  Deborah  Silver  with  the  Center  for  Advanced  Information  Processing,  Rutgers 
University,  is  developing  software  for  repositioning  the  animal  models. 
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3.3  Electrical  Properties 

Determining  permittivity  values  (f'and  cr)  of  various  biological  tissues  is  a  critical 
step  when  calculating  the  SAR.  However,  various  authors'^®’  report  different 
permittivity  values  for  the  same  tissue  type.  This  lack  of  consensus  on  which  permittivity 
values  represent  the  most  accurate  values  of  the  biological  tissues  poses  the  question  of  the 
effect  of  permittivity  value  on  calculated  SAR  values  in  biological  systems.  Since  there  is 
some  variability  in  the  tissue  permittivity  values  of  human  and  animals  and  some 
uncertainty  in  the  measurement  of  the  permittivity  of  tissues,  the  dependence  of  SAR  on 
permittivity  changes  is  important^’. 

The  most  recent  and  comprehensive  permittivity  data  set  includes'”^  more  than  30 
tissue  types  in  the  frequency  range  10  Hz  to  20  GHz  and  on  more  than  20  other  types  from 
10  to  20  GHz.  Three  main  interaction  mechanisms,  each  governed  by  its  own 
kinetics,  determine  the  main  features  of  the  permittivity  spectrum  of  a  biological  tissue. 
Three  main  spectral  regions  known  as  the  a,  p  and  y  dispersions  are  predicted  from  known 
mteraction  mechanisms  (Fig.  10).  These  dispersions  have  been  identified  experimentally  in 
the  frequency  range  from  hertz  to  gigahertz.  The  y  dispersion,  in  the  gigahertz  region,  is 
due  to  the  polarization  of  water  molecules.  The  p  dispersion,  typically  in  the  hundreds  of 
kilohertz  region,  is  due  mainly  to  the  polarization  of  cellular  membranes  that  act  as  barriers 
to  the  flow  of  ions  between  the  infra-  and  extracellular  media.  Other  contributions  to  the  p 
dispersion  come  from  the  polarization  of  protein  and  other  organic  macromolecules.  The 

low  frequency  a  dispersion  is  associated  with  ionic  diffusion  processes  at  the  site  of  the 
cellular  membrane”*. 


«  •  value  (8’)  and  conductivity  (a  in  niS/m)ofthe  muscle  are  shown  as 

a  function  of  frequency  .  Three  major  dispersion  a,  p,  and  y  are  typical  for  all  tissues  and  cells  in 
suspension,  although  magnitudes  and  dispersion  frequencies  vary.  Additionally  smaller  relaxation  effects 
contribute  to  the  high  frequency  -tails  of  a  and  p  dispersions,  as  indicated  by  the  dashed  curves. 


18 


We  use  an  HP  8510  Network  Analyzer  (Fig.  1 1)  to  measure  the  reflection 
coefficient  fi'om  an  open  ended  coaxial  probe  that  is  in  contact  with  the  material  of 
interest  (biological  tissue  or  liquid  media).  These  reflection  coefficient  values  are  then 
used  to  calculate  the  electrical  properties  (permittivity)  of  this  material'*®.  Measurements 
can  be  performed  for  fi-equencies  from  45  MHz  to  26.5  GHz. 


Figure  11.  HP  8510  Network  Analyzer 
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4.0  INTERNATIONAL  EMF  DOSIMETRY  HANDBOOK 


The  5th  Edition  of  the  Radio  Frequency  Dosimetry  Handbook  (RFDH)  will 
expand  on  the  4“’  Edition  by  updating  RF  dosimetry  research  that  has  been  accomplished 
since  1986,  primarily  in  the  fields  of  localized  specific  absorption  rate  (SAR)  and 
temperature  distribution  resulting  fi-om  exposure  to  RF  fields.  The  5*  Edition  will  have 
approximately  two  dozen  chapters  included  under  the  broad  categories  of  general  RF 
principles,  computational  modeling,  analytical/empirical  dosimetry,  experimental 
dosimetry,  and  exposure  assessment.  Dr.  E.  Adair  will  be  writing  the  chapter  on  Thermal 
Responses  to  RF  Exposure,  Prof.  M.  Ziskin  on  Responses  to  RF  Overexposure,  S.Allen 
on  Empirical  RF  Dosimetry  Techmques,  W.  Hurt  on  Basics  of  Electromagnetics,  and  Dr. 
J.  Leonowich  on  Personal  Dosimetry  Indicators.  We  will  contact  about  six  other  experts 
to  write  additional  chapters.  Dr.  P.  Chadwick,  Microwave  Consultants  Limited,  England, 
is  negotiating  vrith  other  experts  to  write  approximately  14  chapters  as  well.  The 
projected  schedule  calls  for  the  authors  to  submit  their  chapters  in  draft  form  starting  in 
2003. 


5.0  SUMMARY 


Knowledge  of  internal  electric  and  magnetic  fields,  induced  current  densities,  and 
SARs  is  of  basic  interest  in  the  assessment  of  biological  effects  and  medical  applications 
of  electromagnetic  fields.  Dosimetry  is  an  important  part  of  any  scientific  effort  to  assess 
the  effects  of  EMF  on  biological  organisms.  Experimental  dosimetry  is  essential  in 
determining  the  internal  fields  and  the  whole  body  or  localized  SAR  values  in 
experiments  with  phantoms  or  animals.  Furthermore,  carefully  performed  experiments 
are  crucial  in  verification  of  theoretical  predictions  and  delineation  of  their  limitations.  In 
the  last  few  years,  considerable  progress  has  been  achieved  in  experimental  and 
numerical  dosimetry.  Today  a  broad  range  of  tools  for  dosimetric  analysis  for  the  wide 
fi-equency  range  is  available.  Because  of  the  heterogeneity  of  the  tissues  and  the 
nonuniformity  of  the  incident  fields,  closed-form  analytical  solutions  are  impossible  and 
computer  methods  are  needed  to  obtain  the  internally  coupled  fields.  In  theoretical 
dosimetry,  FDTD  is  currently  the  most  acceptable  choice  if  digital  anatomical  models  of 
man  and  animals  (highly  inhomogeneous  structures  of  millimeter  resolution)  are  to  be 
analyzed.  The  main  contribution  of  our  research  efforts  to  the  discipline  of  theoretical 
dosimetry  is  identification  of  the  dependence  of  predicted  SAR  in  relation  to  the 
variability  in  permittivity  values  for  different  tissue  types  (muscle,  fat,  skin^  bone 
marrow)  when  using  different  numerical  anatomical  models.  Whole  body  SAR  is  not 
particular  sensitive  to  variations  in  permittivity  values,  while  localized  SAR  values  have 
a  great  dependence  on  them.  The  ratios  in  whole  body  and  localized  SAR  values  are  also 
dependent  on  different  exposure  conditions  (orientation,  firequency). 
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